Luttinger liquid with asymmetric dispersion 
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We present an extension of the Tomonaga-Luttinger model in which left and right-moving particles 
have different Fermi velocities. We derive expressions for one-particle Green's functions, momentum- 
distributions, density of states, charge compressibility and conductivity as functions of both the 
velocity difference e and the strength of the interaction /3. This allows us to identify a novel 
restricted region in the parameter space in which the system keeps the main features of a Luttinger 
liquid but with an unusual behavior of the density of states and the static charge compressibility 
K. In particular k diverges on the boundary of the restricted region, indicating the occurrence of a 
phase transition. 
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In the last years there has been much interest in the 
study of one-dimensional (ID) condensed matter prob- 
lems Specific examples of experimentally realized 
ID structures are: strongly anisotropic organic conduc- 
tors Ql , charge transfer salts |^ , quantum wires Q , edge 
states in a two-dimensional (2D) electron system in the 
fractional quantum Hall (FQH) regime Q and the re- 
cently built Carbon Nanotubes ||]. All these systems 
are no longer described by the usual 3D-like Fermi liquid 
picture. They are believed to belong to a novel, highly 
correlated state of matter known as the Luttinger liquid 
(LL) 01 . Very recently, possible LL behavior in 2D high 
temperature superconductors has also been reported [^. 

From the theoretical point of view the most wide! 
studied ID model is the so-called "g-ology" model 
which is known to display the LL behavior characterized 
by spin-charge separation and by non-universal (interac- 
tion dependent) power-law correlation functions. In par- 
ticular it predicts a momentum distribution function that 
vanishes at pp as n(p) ^ (p — pp)'^'^, where 7 is related to 
the strength of the electron-electron interaction (in the 
free case one has 7 = and n{p) ^ d{pp +p)). One of 
the simplest and yet very useful version of the "g-ology" 
model is the exactly solvable Tomonaga-Luttinger (TL) 
model, which describes left and right-moving electrons 
subjected to forward-scattering interactions | |lO[ |. 

In this Letter we propose a simple modification of the 
TL model in which left and right-moving electrons have 
different Fermi velocities vl and v^. Previous studies 
of LL systems involving more than one Fermi velocity 
are related to an special class of chiral LL and to 
multiband and multichain models [|l2j. Another inter- 
esting problem in which one has different values for vp 
is the interaction between parallel conductors leading to 
the so called Coulomb drag [Q. We want to stress that 
the model we shall study is crucially different from all 
these systems since it is neither a purely chiral LL nor a 
multiband system with symmetric dispersion. Our the- 



ory is formally similar to a recently proposed model for 
the study of spin-orbit coupling in interacting quasi-lD 
systems iQ. These authors, however, concentrated their 
attention on the interplay between velocity asymmetry 
and spin degrees of freedom, whereas here we derive and 
analyze physical consequences connected to the asym- 
metric dispersion only. As we shall see, this point of 
view allows us to obtain some novel non trivial features 
of the system. 

To be specific we start by considering an asymmetric 
dispersion described by the following Hamiltonian 



H = -ih I dx {vR^pl^doolljB. - VLtpLd^IpL) 

+ TrU dx^^j^il:Rijlil:L, (1) 



where i^r^l and -i/'Jj l ^'^^ electron operators and U is 
the strength of the forward-scattering electron-electron 
interaction. In the "g-ology" language we have g2 = t^U 
and (74 = 0. The extension of our results to the general 
case {g4 ^ 0) is straightforward, here we consider this 
particular case in order to keep the discussion as clear as 
possible. We will set ft = 1 from now on. Please note that 
both I'L and vr are positive, and U > corresponds to 
repulsive interactions. This is the case we shall examine 
throughout this work. 

Since the edge states of FQH systems have been suc- 
cessfully described in terms of chiral fermions with drift 
velocities proportional to E/B js) (B is the uniform 
transverse magnetic field and E is an electric field that 
keeps electrons inside the sample ||l^), the model above 
could be experimentally realized by putting together the 
edges of two FQH samples in the presence of different 
fields such that the resulting fractions are also different. 
In such experimental array U represents the strength of 
the interaction between the charge-densities (CD) of each 
fermionic branch. Recent experiments on tunneling be- 
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tween edge states of laterally separated quantum Hall 
effect systems seems to indicate that the experiment 
we propose is indeed feasible. 

One interesting result of this Letter is the appearance, 
due to the velocity asymmetry, of a new available re- 
gion in the space of couplings in which the model (1) 
predicts an anisotropic phase, in the sense that the col- 
lective charge-density modes associated to each branch 
propagate in the same direction. However one has to 
be cautious with this prediction since our computations 
show that (1) is no longer valid in this region. Indeed, as 
we shall see, if we define 



«o = {vr + wl)/2 



and 



VR - VL 
VR + VL ' 



P = U/2vo, 



(2) 



(3) 



we can only trust our model inside the unit circle in the 
e—(3 parameter space (e^-|-/3^ = 1). We then first explore 
the physics described by (1) in the restricted region. We 
shall be specially interested in discussing the cases of con- 
stant asymmetry (e fixed) and constant interbranch in- 
teraction (/? fixed) . In so doing we found a drastic change 
in the behavior of the charge compressibility k where the 
value at zero asymmetry is multiplied by a factor which 
diverges on the transition curve. Studied as a function of 
f3 for fixed e, it first reaches a minimum and then there is 
a strong enhancement as /3 ^ /3o = ^/vrVl/vq, in oppo- 
sition to the monotonous decay present in the e = case. 
A similar change of behavior is present in the density of 
states (DOS) function. For /3 and e sufficiently small one 
recovers an ordinary LL system (wq ^ vf,v- ^v+). 

Outside the restricted region there is a change in the 
sign of one of the "plasmon" velocities, accompanied by 
a dramatic change in the behavior of the Green function 
for that branch, which now diverges at long distances. 
This, together with the fact that k becomes negative, 
bring out that the model suffers some kind of instability 
in this region. It is important to stress that this region 
is absent tor vr — v^. 

We have studied the system by using functional 
bosonization techniques p^ . This amounts to defining 
fermionic field operators in the Heisenberg picture. We 
then have a field-theoretical, Lagrangian formulation of 
the model. This, in turn, allowed us to obtain an action 
describing the dynamics of the bosonic collective excita- 
tions of the system. Using this action one can easily com- 
pute the dispersion relations for the CD oscillations. For 
short-range, constant electron-electron potentials, these 
dispersions are linear, with velocities given by v± — wo?y±j 
and 



From this equation one sees that the propagation of the 
collective modes takes place for /3 < 1. It becomes appar- 
ent that, in contrast to the usual answer for a TL model 
with vr — V]^ and 174 = 0, here one has two different 
velocities v+ and V- for the propagation of left and right 
CD modes. Moreover, one of the velocities w+ or v- goes 
to zero as the interaction and the asymmetry approach 
the curve + f]'^ = 1 and changes its sign beyond that 
curve, as anticipated above. If one keeps (3 fixed this 
change of sign occurs for e — €q = \J\ — fS^ . At this 
point, as we will see, there is a divergence in the charge 
compressibility and in the DOS, which suggests that a 
phase transition takes place. 

Now, in order to get an insight into the physical 
consequences of the velocity difference, we compute 
single-particle quantities: the Green function Gr(x,t) — 
(V.(x,t)V'J(0,0))^^^^ with r = +(i?),-(L), the mo- 
mentum distribution function, the spectral function 
Pr{q,i^) given by 



Pr('Z,^) = --ImG(-*)(g,a;), 



and the DOS defined as 



(5) 



(6) 



In these equations, G^™'' (q, lo) is the Fourier transform 
of the retarded Green function: 



Gf-^"') {x,t)^ lOit) {{M^,t),i,l (0, 0) } ) . (7) 

For the normal phase -|-e^ < 1), the Green function 
at T = is given by 



Gr{x,t) 



27ra \q; — ir{x — Vrt) 



7+1 



a -f ir{x — v^rt) 



(8) 



where a is an ultraviolet cutoff. The constant 7 has the 
usual expression in terms of the stiffness constant K, 7 = 
{K+K^^ — 2)/4, but in this constant vp must be replaced 
by the average fo: 



K 



vo-U l-(3 



(9) 



vo + U ]l 1 + P 
In the outer region + > 1), tor vr > vl, we get 

1 [a ~ i{x — v^t)]^ 



(4) 



GR{x,t) 



27r [a - i{x - v+t)] 



7+1 



(10) 
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and 



— 1 [a — i{x — v_t)] 
^■Ko? \oL — i{x — v+t)] 



7 + 1 



(11) 



From these results one obtains a momentum distribu- 
tion of the Fermi type for the right branch, i.e., nfi{p) ^ 
9{p — pfi). However, the situation is very different for ul. 
Indeed, when taking the appropriate t — > 0~ hmit in or- 
der to employ the usual definition of one finds that 
the correlator increases linearly with distance, instead of 
having the decay (typical of 3D-like systems) as is 
the case for the right branch, or the x"'^^"*"^-* behavior 
that yields the LL result in the "normal" LL region. This 
leads to til ~ 5'{p + pl)- The appearance of this diver- 
gent (at long distances) left correlator, together with a 
momentum distribution which is not a positive definite 
quantity, are clear indications that the model given by 
(1) is unphysical beyond + = 1. (For vl > vr the 
corresponding left and right behaviors are exchanged). 
We then conclude that the model given by (1) yields sen- 
sible results for + 0^ < 1 and from now on we will 
restrict our study to that region. 

From equation (|^) one can calculate the spectral func- 
tion pr{q,uj). We obtain, as in the symmetric case, only 
one singularity in the positive frequency sector and one 
in the negative sector, as expected for spinless systems. 
The function diverges at those points as 



X [9{uj — v^q)9{uj — v+q) + 9{v^q — uj)9{v+q — to)]. 

(12) 

The exponents do not depend on e whereas the position 
of the singularities does. 

Concerning the DOS we get 



N{uj) = No{lu) 1 



-1-7 



(13) 



with 



Noitu) 



(w/wo) 



27 



7rvor(27 + l)(l-/?2)7+i/2^ 



(14) 



where ujo = vo/a and F is the Gamma function. We 
see that as the asymmetry is increased, the DOS grows 
from its value at e = 0, and diverges at the point e — eo. 
We want to stress that in systems with several spectral 
branches as multicomponent Tomonaga-Luttinger model 
and spin-polarized Luttinger liquids a growing 
of the DOS as increasing the velocity difference between 
spectral branches is also observed. This disagrees with 



the result obtained in the system with spin-orbit coupling 
0-. 

Using standard linear response theory one can express 
the conductivity as an integral of the retarded current 
density correlation function. At this point one has to 
recall that the naive definition of the current j = vrpr — 
vlPl does not satisfy the continuity equation ||l^. This 
choice for j leads to a frequency-dependent conductivity 
that diverges on the unit circle + (3'^ = 1 and becomes 
negative for + P"^ > 1. However, as shown in ||T^ it is 
indeed possible to build a physical current starting from 
the continuity equation. Extending this procedure for 
the present e 7^ case we obtain jptys = (1 — f3){pR — 
Pl) + vo^{pr + Pl)- Using this expression we were able 
to get the frequency-dependent conductivity a as: 



cr(w) 



^^0 (1 - P) 



11 iO 



(15) 



which is independent of e. 

Let us now consider the static charge compressibility of 
this system, defined as k = {{pr + pl){pr + Pl)) (q,^) 
for cj ^ 0. A straightforward computation yields 



Kq 



with 



Kq 



TTVO (1 +^)' 



(16) 



(17) 



where one sees the divergence that takes place, at fixed 
/3, for e — eq. This is similar to the behavior of the DOS, 
although both functions diverge with different exponents. 
Note that beyond cq one obtains negative values for the 
compressibility, a further indication that the model is not 
valid in that region. In Figure 1 we show the dependence 
of K on e for different values of f3. We see that the asym- 
metry enhances the compressibility. Of course, it is also 
possible to study k as function of the coupling /3, for 
a fixed asymmetry. This is depicted in Figure 2 where 
one sees that in drastic departure from the symmetric 
case, which displays decreasing k for increasing /3, now 
K reaches a minimum and then grows without bound as 

The critical behaviour of the static charge compress- 
ibility and the DOS at eo together with the "freezing" of 
one of the spectral branches (w+ 0) is an indication 
that a phase transition involving CD degrees of freedom 
takes place on the boundary 0^ -\- = 1. A similar tran- 
sition related to spin variables was also found in . 

In summary, we have presented a simple modification 
of the usual TL model, in which left and right-moving 
particles have different Fermi velocities. By using func- 
tional bosonization methods we computed the dispersion 
relations of the underlying bosonic collective modes of the 
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FIG. 1: Static charge compressibility as function of e, for 
f3 = 0.1,0.6,0.9. 
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FIG. 2: : Static charge compressibihty as function of /3, for 



system. We showed that the velocity asymmetry gives 
rise to some remarkable features. We have found that 
the values for the electron-electron coupling /3 and the 
asymmetry e are restricted to lie inside the circumfer- 
ence + = 1 in order to guarantee the stability of the 
LL described by (1). In this region, the DOS N{uj) and 
the static charge compressibility k display a big devia- 
tion from the standard LL behavior. For a fixed asym- 
metry one sees that now k reaches a minimum and then 
grows without bound as /3 ^ /Jq- On the boundary of 
the restricted region the velocity of one spectral branch 
goes to zero. Also k and DOS diverge when approach- 
ing this critical boundary, indicating the appearance of a 
new phase. 

To conclude we would like to stress that the idealized 
model we present here could be approximately realized 
by allowing the interaction between the edge states of two 



FQH plates. Since the drift velocities of the correspond- 
ing chiral LL's are proportional to E/B, one could have 
vb. 7^ vl by conveniently tuning up the corresponding 
values of these fractions. 
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